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\ It is well known that the Onsager reciprocal relations give important identities between 

[t I , transport coefficients for heat and charge conduction in the linear response regime?"^ For 

CN ■ example, the Peltier coefficient 11 is related to the Seebeck coefficient S" in a simple form as 
S = n/r, where T is the temperature. The Onsager-Casimir symmetry, which is the extension 

\ of the Onsager relation in the presence of an external magnetic field, is also derived from 

^ , the microscopic reversibility.'^ Recently, a general relation called the steady-state fluctuation 

^ ■ theorem has been derived by the same concept of the microscopic reversibility.'^'^ It gives an 



O 



identity equation on the probability of the entropy production AS" during time r as 

P(AS) 



In 



P(-AS) 



A5, (1) 



\ in the asymptotic limit t ^ oo. While this expression reproduces the ordinal Onsager relations 

O , in the linear response regime, it provides useful information even in the far-from-equilibrium 
regime. 

CN ■ The importance of the Onsager-Casimir relation on coherent quantum transport of meso- 

l/^ \ scopic devices is known for long time.'^'^ Recently, Saito and Utsumi have derived a quantum 
00 \ 

, version of the fluctuation theorem, and applied it for derivation of universal relations among 



nonlinear transport coefficientdSJ by means of the full-counting statistics in the Keldysh for- 



p 

■ malism.l^^^ In Ref. [8l however, its physical meaning on thermoelectric effects has not been 



addressed in detail. 

In this note, we discuss extension of the Onsager relation between the Peltier effect and 
^ . the Seebeck effect to nonequilibrium steady states in mesoscopic devices. For simplicity, we 

■ consider a two-terminal setup, though extension to a multi-terminal setup is straightforward. 

We consider the Hamiltonian H = Ylir=h R + Hd + Ht, where 

Hr ='^^erkali.ark, (2) 



m 



Hd=J2 ^i4di + J2 + h.c.) + Hi, (3) 

i=l (iJ) 
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Ht=Y^ trki{d\ark + h.c). (4) 

rki 

Each of two leads is described by {r = L,R), where ark is an annihilation operator of 
an electron with a wave vector k. A mesoscopic device is modeled in a general form by 
Hd consisting of an arbitrary number of local energy levels labeled by (1 < i < m), 
where di is an annihilation operator at the ith site. In this model, electron hopping tij and 
electron-electron interaction (denoted by Hj) between arbitrary pairs of sites are assumed-l^ 
A coupling between the leads and the mesoscopic device is described by Hj- with electron 
hopping trki- An external magnetic field B is introduced by the Peierls phase on the hopping 
elements as t^-ki = \trki\ exp(i(^j.fci) and tij = \tij \ exp{i(f)ij) where cprki and cpij are odd functions 
of the magnetic field: (j){—B) = —(j){B). In this note, we use a unit h = ks = e = 1. 

We introduce a cumulant generating function (CGF) T{xc,Xe', B) for the steady state, 
where Xc and Xe are counting fields for charge and heat current. Current operators are defined 
as Ic = Nl = ^Wl, Ht] and Ig = -^L = '^[Hl, Ht]^ where Nl is a number of particle of the 
left lead. The CGF generates cumulants by the derivatives with respect to the counting fields. 
For example, the charge current and noise are generated as 



m = {Ic) 



diXc 



(5) 

Xc=Xe=0 



ill)) = {ll)-{Ic?- 



(6) 

Xc=Xe=0 



{dix. 

respectively. General relations among nonlinear transport coefficients are derived from a sym- 
metry relation 

/•(Xc, Xe; B) = T{-Xc + iAc, -Xe + iAe] -B), (7) 

which is a consequence of the microscopic reversibility.'^ Here, Ac = PifJ-L — I^rI^r and 
Ae = — — Pr) are affinities. All the cumulants are expanded with respect to Ac and Ae as 



oo oo rki,k2 



(9) 

Ac=Ae=0 



{{I'c'il')) = Y.Y.TnV^'cA[\ (8) 

'^''^ dA['dAi^ 
By symmetrizing and anti-symmetrizing the transport coefficients with respect to B as 

± (B) = Lit (B) ± Llf^ i-B) , (10) 
general relations are derived from Eq. d?]) a^ 

V ( 1 \ni+n2+ki+k2 t ki+ni ,k2+n2 \ 
^ ^~^) ^h-ni,l2-n2,±- ^^^) 
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We can show immediately that Eq. (jlip includes the Onsager relation in the linear response 
region as follows. First, L'l^^f^^{B) are related to ordinary observables by changing affinities to 
the external-bias variables, Ac and Ae to the bias voltage V = fii — t^R and the temperature 
difference AT = Tl — Tr. (For example, the linear conductance of charge current, the charge 
current noise at equilibrium and the Seebeck coefficient are written as G = l|'q/T, Sjj = Lqq 
and S = Lqi/TL\'q, respectively.) We note that the linear response of charge and heat 
currents are expressed by Ic = lI'°q{B){V/T) + lI'°^{B){AT/T^) and 1^ = L°{^q{B){V/T) + 
Lq'\{B){AT/T'^)^ Then, specific equations between coefficients satisfying N = ki + k2 + li + 
^2 = 2 lead to the Onsager relation L^'^{B) = l1'^{—B). 

Next, we prove that the magnitude of the nonlinear Peltier effect can be determined only 
by information of charge current measurements without heat current measurement; we show 
that all the coefficients for heat current appearing in the expansion 

oo ^0,1 



under the isothermal condition {Ae = 0) can always be rewritten by the transport coefficients 
of the higher-order charge curreni 
general equation p!T|) . we obtain 



of the higher-order charge current cumulant L'^^f. By substituting k2 = and I2 = ^ into the 



N 

where M is a matrix whose matrix elements are given as 



-ly (0<i<j<N) 
M,, = <( \J-^J' y - -J- , _ ^^^^ 

(0 < i < i < iv) 

By utilizing MijMjk = 6ik, we obtain 

N 

^n]o,± = T E ^ofc-^^^°_fc,i,±- (15) 

k=l 

We notice that Eq. (|15|) relates the nonlinear response of the heat current under Ae = to 
that of the charge current, current noise, and higher cumulants up to the linear response in 
Ae- This means that the magnitude of nonlinear Peltier effect can be evaluated without direct 
measurement of the heat current. This result can be regarded as a nonequilibrium extension 
of the Kelvin-Onsager relation. 

Finally, we discuss the lowest-order nonlinear correction of the Peltier coefficient defined 

by 

uiAc) = Y = n^°^ + n(i)/c + o(/2). (16) 

under the isothermal condition {Ae = 0). The linear-response term satisfies the Kelvin- 
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Onsager relation 11^'^) = TS^^\—B) as already mentioned, where S^'^\B) is the ordinal 
Seebeck coefficient defined in the linear-response regime. The first-order correction 11^^^ is 
formally expressed by using the definition of L^^f^ as 



n(i) 



2L\l{BY \L%{B) L\l[B)l ' 
By utilizing the general relation in Eq. (|lip for N = 2, the correction is rewritten as 

n(i) _ S{-B) ( L'^liB) 



T 2 \l\1{BY 

lZ{-B)-2L\1{-B) \ 
L\l{B)L]il{-B) )■ 

Thus, though the expression is a little complicated, one can find that the nonlinear correction 
of the Peltier coefficient can be calculated only by coefficients L'^^f^, which needs no heat- 
current measurements. In a similar way, higher-order correction of the Peltier coefficient (e.g. 
n'^^-') can be expressed by higher-order cumulants of the charge current (e.g. skewness). 

We proved that the magnitude of the nonlinear Peltier effect under the isothermal con- 
dition can be evaluated from transport coefficients for charge current cumulants (current, 
current noise, and higher cumulants) up to the linear response to a thermal bias. Inversely, if 
precise measurement of both charge and heat currents is possible, thermoelectric effects in the 
nonlinear regime can be used for an experimental proof of the fluctuation theorem in quan- 
tum systems, as performed in a recent experiment by means of simultaneous measurement 
of conductance and shot noise.'S By combining recent progress in measurement techniques 
of thermoelectric effects,'^ the present result will provide an important key to understand 
nonequilibrium thermoelectric effects in mesoscpic systems. 
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